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^ , A unified framework for coupled Navier-Stokes / Cahn-Hilliard equations is 

h : 

developed using, as a basis, a balance law for microforces in conjunction with 
constitutive equations consistent with a mechanical version of the second law. 
As a numerical application of the theory, we consider the kinetics of coarsening 
for a binary fluid in two space-dimensions. 
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I. INTRODUCTION 



The Cahn-Hilliard equation [T|] [] 



if - = mA [f'(<p) - aAip] 



(1) 



is central to materials science, as it characterizes important qualitative features of two-phase 
systems. This equation is based on a free energy 



with f((p) a double- well potential whose wells define the phases, and leads to an interfacial 
layer within which the density ip suffers large variations. 

The standard derivation of the Cahn-Hilliard equation begins with the mass balance 



and the constitutive equation 



1 Notation. Tensors are linear transformations of i? 3 into i? 3 and are denoted by upper-case 
boldface letters. Vectors may be viewed as 3x1 column vectors and tensors as 3x3 matrices. 1 
denotes the unit tensor; a ® b, the tensor product of vectors a and b, is the tensor defined by 
(a (g) b)u = (b • u)a for all vectors u; A T is the transpose of a tensor A; tr A is the trace of A; the 
inner product of tensors A and B is defined by A • B = tr(A T B). 

The divergence, gradient, and Laplacian of a field <p = </?(x, t) are denoted by grady?, divcp, and 
Atp. For a vector field u(x), gradu(x) is the tensor with components dui/dxj (i = row index, 
j = column index). The divergence of a tensor field A(x) is the vector field with components 
J2j9Aij/dxj (i=row index). 

The derivative of a function / of a scalar variable (not time) is denoted by a prime: /'. The partial 
derivative of a function $(a, b,c,...,d) (of n scalar, vector, or tensor variables) with respect to b, 
say, is written db&(a, b,c, . . . ,d). 




(2) 



— div h 



(3) 
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h = — m grad/Li, (4) 

which relates the mass flux h to the chemical potential \i. The presence of density gradients 
renders (0) incompatible with the classical definition of \i as the partial derivative of ip with 
respect to ip; instead \i is defined as the variational derivative /i = 8^ /Sip of the total free 
energy 

ty((p) = / ip(ip, gr&d(p)dv, B = underlying region of space. (5) 

Since 

6*/5<p = f'(<p) - aA v , (6) 

this yields the Cahn-Hilliard equation. 

The major advances in nonlinear continuum mechanics over the past thirty years are 
based on the separation of balance laws (such as those for mass and force), which are 
general and hold for large classes of materials, from constitutive equations (such as those 
for elastic solids and viscous fluids), which delineate specific classes of material behavior. In 
the derivation presented above there is no such separation, and it is not clear whether there 
is an underlying balance law that can form a basis for more general theories. 

An alternative derivation of the Cahn-Hilliard equation [[J is based on a balance law 
for microforces ||, defined operationally as forces whose working accompanies changes in 
if. These forces are described by a (vector) stress which characterizes forces transmitted 
across surfaces, and a (scalar) body force ir, which represents internal forces distributed over 
the volume of the material. The basic hypothesis is the microforce balance 

/ £ ■ nda + [ irdv = (7) 
JdR Jr 

for each control volume R, with n the outward unit normal to dR. 

Here we study the isothermal motion of an incompressible binary fluid, with one of the 
constituent densities serving as the order parameter. A chief assumption is that the relative 
momenta and kinetic energies of the constituents are negligible when computed relative to 
the gross motion of the fluid. This allows us to base the theory on: 
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• balance of mass for the order parameter; 

• the microforce balance (|7]); 

• balance of momentum for the macroscopic motion of the fluid; 

• a version of the second law in the form of a global free-energy inequality in which the 
microscopic working 



joins the working of the macroscopic stresses; and 
• constitutive equations presumed compatible with this version of the second law. 

In standard theories of diffusion the chemical potential is given, constitutively, as a 
function of the density, but here we consider systems sufficiently far from equilibrium that 
a relation of this type is no longer valid; instead we allow the chemical potential and its 
gradient to join the stretching tensor, the density, and the density gradient as independent 
constitutive variables. Interestingly, the second law reduces the free energy to a function 
ip = ip(ip, grad<£>) of at most tp and grady, where gradp here denotes the spatial (Eulerian) 



A chief result of our theory is a constitutive equation giving the macroscopic (Cauchy) 
stress as a classical Newtonian stress plus a term 




(8) 



gradient. 



a gr&dp ® grady?, 



(9) 



which represents capillarity. Our final results consist of the coupled Navier-Stokes/Cahn 



Hilliard equations 



p [vt + (gradv)v] 



gradp + z/Av — ct(A(/?)grad<yj, divv = 0, 




(10) 



2 A term of this type appears in Korteweg's theory of gradient fluids (cf. Truesdell & Noll 
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We also develop equations for the motion of immiscible fluids. Here our use of an 
order parameter endows the interface with capillarity and yields a theory in which the free- 
boundary need not be explicitly tracked; as this theory does not allow for a phase change, 
the interface is transported with the material. 

The model defined by (|T0D, which is referred to as "Model H" in the literature on crit- 
ical phenomena ||, was used by Siggia, Halperin, & Hohenberg || to study behavior at 
the critical points of single and binary fluids. Later extensions of the model were used to 
analyze the superfluid transition of 4 He |7j as well as the decay of critical fluctuations in 
simple fluids |§ and polymer solutions |J under an imposed shear flow. A chief purpose of 
these studies was to analyze the effect of reversible modes (i.e., those corresponding to the 
coupling terms (grady?) • v and a(Ay?)grady?) on the dynamical critical behavior of otherwise 
purely dissipative models, with emphasis placed on the behavior of hydrodynamic transport 
coefficients near the critical point. More recently, (|10"D was used to study other nonequilib- 
rium phenomena below the critical point; for example, Kawasaki & Ohta |10[ and Koga & 
Kawasaki [|TT|] use JIHD to study hydrodynamic effects during spinodal decomposition. 

The applicability of fllUp to studies of the type discussed above is generally based, not 
on first principles, but instead on the observation that a convective term (grad<y?) • v in 
the Cahn-Hilliard equation should be accompanied by a term a(A<^)gradv? in the Navier- 
Stokes equation to ensure that the asymptotic probability-distribution functions obtained 
as t — > oo coincide with those given by standard equilibrium calculations. In view of 
the phenomenological nature of this formalism, it seems natural to ask whether there is a 
derivation of ([10]) within the framework of nonlinear continuum mechanics. Here we provide 
such a derivation. 

In addition, we present computational results appropriate to a range of parameters in 
which sharp interfaces (large, localized changes in <p) exist. 

From a computational standpoint our approach is similar to VOF (volume-of-fluid) meth- 
ods |I2J. There a color function, which plays a role similar to that of if, is passively advected 
by the flow. More recently the VOF method has been extended to include capillary effects 
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T3| , |T4|| leading to a modified Navier-Stokes equation equivalent to the first of (|10D. These 



models, however, do not endow the color function with any physical significance, and dissi- 



pative relaxation as in the equation for (p in (10) is absent. In comparison, not only does the 



model based on ( |T0| ) include capillary effects in the modified Navier-Stokes equation, it also 
explicitly includes other dissipative processes up to the length scale of the smeared interface; 
it can therefore describe events such as break-up and coalescence, with the proviso that the 
physical behavior at those length scales is determined by the choice of free energy ijj. 

II. BASIC LAWS 

We consider the isothermal motion of an incompressible binary fluid whose total density 
is constant. A basic assumption of our theory is that the momenta and kinetic energies of 
the constituents are negligible when computed relative to the gross motion of the fluid. 

A. Kinematics 

We write grad and div the spatial gradient and spatial divergence. Incompressibility 
then requires that v, the gross velocity of the fluid, satisfy 

div v = tr (gradv) = 0, (11) 

a constraint that renders the stretching, 

D = - (gradv + gradv r ) (12) 

traceless. 

We will use both the spatial time derivative (...)( and the material time derivative (•••)'; 
these are related through 

(...)• =(...)* + [grad(...)]v. (13) 
We will also use the following identity for scalar fields (p, 

grad(<£>') = (grad<£>) ' + (gradv) T grad</9. (14) 
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B. Balance of mass 



{/»•*} 






JR 



We write pi, for the density and for the mass flux of constituent b (= 1,2), with 
measured relative to the gross motion of the fluid. Then 

Pi + P2 = P (15) 

with p, the total density of the fluid, constant; consistent with this constraint, we assume 
that 

hi + h 2 = 0. (16) 

Throughout the paper R denotes an arbitrary control volume (fixed region of space), 
with n the outward unit normal to dR. Given a field $, 

$ m dv, (17) 

where the dot signifies the time derivative following the material currently in R, and where 
the integral and time derivative commute by virtue of balance of mass, since the total density 
of the fluid is constant. 

Balance of mass requires that 

jy p&G?i>j =—J hft • n da (18) 

for each constituent b. Because of fll5D and fllBD , one of the two relations (|18D is redundant; 
we therefore let 

<P = Pi h = hi, (19) 
and restrict attention to the local balance 

<p m = -divh (20) 

for the first constituent. 
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C. Balance of momentum 



We write T for the stress tensor associated with the macroscopic motion of the fluid. 
Then, since we neglect the relative momenta of the constituents, the balance laws for linear 
and angular momentum have the standard form 

f Tnda = I f p\dv\ f x x Tnda = j / (x x p\)dv j (21) 

J dR y.JR J J dR y.JR J 

for each R, or equivalently, 

divT = pV T = T T (22) 
It is convenient to introduce the extra stress S and the pressure p defined by 

S = T + pl p = --(trT)l. (23) 
As a consequence of incompressibility, p is indeterminate. 

D. Order parameter. Microforce balance 

We assume that the microscopic behavior of the fluid, as manifested in the diffusion of 
its constituents, is described by a scalar order-parameter u and concomitant microforces 
whose working accompanies changes in u. Since working is characterized by terms of the 
form { microforce times W }, microforces are here scalar quantities. Precisely, we assume 
that the microforces are described by a (vector) stress £ -whose traction £ • n characterizes 
forces transmitted across oriented surfaces of unit normal n- in conjunction with a (scalar) 
body force tt, which represents internal forces distributed over the volume of the fluid; these 
forces are presumed consistent with the microforce balance [H 

/ i-nda+ [ irdv = (24) 

JdR JR 

for each control volume R, or equivalently 

div£ + tt = 0. (25) 
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A basic hypothesis of the theory is an identification of the order parameter uo with the 
constituent density ip: 

Lu = cp. (26) 



E. Second law in the form of a dissipation inequality 

We restrict attention to isothermal behavior and therefore consider a mechanical version 
of the second law: for each control volume R, the rate at which the energy of R increases 
cannot exceed the working on R plus the rate at which energy is transported to R by 
diffusion. Let /x& denote the chemical potential of constituent b, and let 

li = iii - /i 2 (27) 

Then 

— ^2 J /^fehfe ' n <^ a = ~ J A*h ' n da (28) 

6=1,2 ®^ ®^ 

represents energy carried into R across dR by diffusion, while 

/ T n- vda, [ u>' £ • nda (29) 

JdR JdR 

gives the working of the macroscopic and microscopic stresses. Thus, since we neglect the 
relative kinetic energy of the constituents, the appropriate form of the second law is the 
dissipation inequality 

\ [ (tp + k)dv\ <[ Tn-vda+f if' £-nda- [ uh-nda, (30) 

Ur ) JdR JdR JdR 

where k = pv 2 /2. The body force ir does not contribute, since it acts internally to R. 
By (0), fll|), (HD, and (H), © has the local form 

i/j ' - T • gradv + (vr - /i)ip ' + £ ■ grad(<^ ' ) + h • grad/x < 0. (31) 

Combining (|3lD with (|l^), (|22|) , and (|23|) , we are led to a local dissipation inequality 



ip' — [S + (grady?) (g> £] • gradv + (tc — fj,)ip ' — £ ■ (grady?) ' + h • grad/i < (32) 

that will form a basis for our discussion of constitutive equations. 

The negative of the left side of (|32|) represents the dissipation V, as its integral over R 
is the right side of (|30|) minus the left. Thus, for motion in a container B with both v = 
and if £ • n = on dB, 

lj(ip + k)dv\ =-Jvdv<0 if h-n = on dB 
jy (ip — /lop + k)dv\ = — J Vdv < if \x = /io = const, on dB. (33) 

The thermodynamic development therefore yields natural Lyapunov functions for certain 
classes of flows. 



III. CONSTITUTIVE EQUATIONS. RESTRICTIONS IMPOSED BY THE 

SECOND LAW 

In standard theories of diffusion the chemical potential is given, constitutively, as a 
function of the density, but here we consider systems sufficiently far from equilibrium that 
a relation of this type is no longer valid; instead we allow the chemical potential and its 
gradient to join the stretching, the density, and the density gradient in the list of constitutive 
variables. Precisely, we consider constitutive equations of the form 

if; = $(D, (p, grady?, /i, grad/i), 
S = S(D, if, grady?, //, grad/i), 

I = £(D, <p, grady?, /i, grad/i), (34) 
7r = 7r(D, (p, grady?, /x, grad/i), 
h = h(D, yj, grady?, /i, grad/i), 

with each of the constitutive functions isotropic, since the material is a fluid. To avoid 
notation such as ^^^(D, ip, grady?, //, grad/i) for the partial derivative with respect to 
grady?, we write 
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b = grady?, s = grad/i. 



(35) 



Not all constitutive relations of the form (34) are admissible, as without further restrictions 
(|34l) will violate the dissipation inequality (p2|) . To determine the requisite restrictions we 
choose arbitrary fields v, ip and p, for the velocity, order parameter, and chemical potential, 
and use (|34] ) to compute a constitutive process consisting of v, (p, p and the fields if}, S, ir 
and h; such a constitutive process will satisfy (|32|) if and only if 



<9dV(- 



D' + 



d^{... )+%{...) -p <p- + <h,c(...) £(...) 



■ b' - 



S(. . .) + b ® £ • gradv + 



d^{...) /i- + a^(...) -s- +h(...)-s<0 



(36) 



where we have written (...) as shorthand for the list of independent constitutive variables: 



vD,v?, gradv?, /i, grad/i). 



(37) 



It is possible to find fields v(x, t), ip(x, t), and /i(x, t) such that v, gradv, (gradv) ' ,ip, ip' , b = 
grad^, b ' = (gradp) • , p, p ' , s = grad/i, and s ' = (grad/i) ' have arbitrarily prescribed 
values at some chosen point and time.| Thus, since the terms D • , ip ' , b " , // " and s ' appear 



linearly in (|36|) , it follows that drt^jj = 0, d^ip — fj, — tt, d^ijj = 0, d s tp = 0, and d^ip = We 
are therefore led to the following constitutive restrictions: 



It is tacit that there are external mass supplies and forces available to ensure satisfaction of 
the balance laws for arbitrary choices of the fields v, cp and [i. We chose not to introduce external 
fields (source terms) since they tend to complicate the discussions and since it is only here that 
they are required. Specifically, we need an external mass supply r and external body forces b and 
7 such that @, <^), and (|2|) become ip' = -divh + r, divT + b = p\' , and div£ + vr + 7 = 0, 
and such that the dissipation inequality ( |30[ ) contains the additional term f R {b • v + ip ' 7 + pr} dv 
on the right hand side. Allowing an external supply for each balance law is an assumption now 
standard in continuum mechanics; such assumptions are, of course, tacit in derivations dependent 
upon arbitrary variations of a field. 
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(i) the free energy and microstress are independent of D,/i, and grad/i, and are related 
through 

£((p, grad(^) = d h ijj(ip, gradp); (38) 



(ii) the internal microforce is independent of D and grad/z and represents a nonequilibrium 
contribution to the chemical potential: 



vr(v9, grady?, fi) = n - d v ijj(ip, gr&dtp); 



(39) 



(iii) the following inequality must be satisfied for all values of the arguments: 



h(. . .) • grad/i — S(. . .) + b ® £(</?, grady?) • gradv < 0. 



(40) 



Since the skew part of gradv can be chosen arbitrarily and independently of D, and since 
S is symmetric, £ = £(ip, b) satisfies 



b®£ = £®b, (b-^)b = |b| 2 ^ 



(41) 



and is hence parallel to b. Further, the isotropy of £ = £(<p, b) implies that £ = when 
b = 0; hence there is a scalar function a(ip, b) such that 



(42) 



Let 



P = P(...) = S(...) + a(^,b) 
then, since trD = 0, (|40|) reduces to 



bob — -|b| 2 l 

3 1 1 . 



(43) 



• s — P(. . .) • D < 0. 



(44) 



P in fl4"4|), represents a thermodynamic stress, as it is conjugate to the stretching D. 

Any set of constitutive equations consistent with the restrictions (ff8D, Q3"9"D, and 
([14]) will be consistent with the dissipation inequality (|32|) . Our purpose here is not to 



12 



develop the most general theory possible, but rather to develop a theory that couples the 
essential features of the Cahn-Hilliard and Navier-Stokes equations. Newtonian fluids have 
stress linear in D, while Cahn-Hilliard diffusion has mass flux linear in grad/i and free 
energy quadratic in grady?. Guided by these theories and by ([44]), we now assume that the 
thermodynamic stress P, and the mass flux h have the specific forms 

P = 2v((p)D, h = -m(y?)grad/i, (45) 

with mobility m(ip) and viscosity z/(y?) nonnegative, and, appealing to d42|) , that the free 
energy ip has the form 

^ = /(^) + i«(^)|grad^| 2 (46) 

with a(ip) nonnegative. Then 

T = —pi + 2u((p)H — a(y?)grady9 ® grady?, 

$, = a((p)grad(f, (47) 
7T = H - f((p) - ^a'(<^)|grad^| 2 , 

where, for convenience, the pressure p has been replaced by p+ (a/3)|grad<^| 2 . If we assume, 
in addition, that a, m, and v are constants, then 





f{ip) + ^a|grady?| 2 




T = 


—pi + 2z/D — a grady? ® grady? 




t = 


a grady?, 


(48) 


71 = 


M - f\v) 




h = 


—m grad/i. 


(49) 



The term a grady? ® grady? in the stress gives rise to normal stresses in the absence of flow; 
this term should represent surface tension and could effect boundary conditions at a free 
surface in a nonstandard manner. 
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Further, and what is most important, the microbalance ( p5|) and the constitutive relation 
for 7r yield Cahn's formula for the chemical potential 

f A = f'( ( p)-aA<p. (50) 

The basic PDEs of the theory follow upon substituting (pE8|) and flSDD into (pD[) and (p2|); 
the results, after replacing the original pressure p by p— (a/6)|grady?| 2 , consist of generalized 
Navier-Stokes equations 

pv' = — gradp + i/Av — a(A<p)grady?, divv = 0, (51) 

coupled to the Cahn-Hilliard equation 

(p- -mA[/V)-aA4 (52) 

The latter equation depends explicitly on the flow velocity, since the material time derivatives 
have spatial forms: 

<£' =<Pt+ (grady?) • v, v* = v t + (gradv)v. (53) 

An interesting alternative form for these PDEs mentions the chemical potential explicitly: 

pv" = — gradp + z/Av + fi grady?, divv = 0, 
(p m = mA/i (54) 

V = f'(<p) - aAl P> 

where the original pressure p has now been replaced by p — (ct/6)|grad<p| 2 + f((p). 

IV. IMMISCIBLE FLUIDS DESCRIBED BY AN ORDER PARAMETER 

With but minor modifications the theory presented above can be applied to immiscible 
fluids. Here the order parameter <p has a given value in each fluid, with the fluid interface 
defined by the variation in if between values. Since the fluids are immiscible, the interface 
should be transported with the fluid; we therefore require that 
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<p- = 0, 



(55) 



so that ip is constant on streamlines. Further, we let p = p((p), allowing the fluids to have 
different densities; this dependence on p is consistent with incompressibility, since, by (p5|), 

piipy =o. 

The microforce balance is no longer relevant; the basic equations are the classical balances 



(plf) and the dissipation inequality 



Up + k)dv> < Tn • vda k = -pUpW 
r J JdR 2 



(56) 



and this leads to the local dissipation inequality 



ip ' — S • gradv < 0. 



As constitutive equations we assume that 



(57) 



-0 = ^(D,^,grady>), S = S(D, <p, grad^), 



(58) 



and proceeding as before, consistency with ([57]) leads, by virtue of (|I4| ) and (|55|), to the 
inequality 



<9 D ^(D,v?, b) 



D 



S(D,v9,b) + b® d h ip(D,p,b) • gradv < 



(59) 



rendering the free energy independent of D. Finally, restricting attention to energies that 
are independent of <p and quadratic in grad<£>, we are led to the constitutive equations 



i> = -ct|grad<^| 2 , 



-pi + 2v(ip)~D — a jgrady <8> gr&dp — -|grady9| 2 l| . 



(60) 



(A dependence of ip on cp does not alter the resulting PDEs.) 

The basic equations, after replacing p by p — (a/6)\gra.dp\ 2 , therefore have the form 

p(<p)v' = — gradp + v((p)Av + 2z/(<y?)Dgrad<£> — a(Ap)gradp, 

divv = 0, tp' = 0. (61) 
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These equations may be useful in the study of the interface between immiscible fluids; the 
order parameter automatically tracks the interface, here a layer, and provides a capillarity 
term a(A<£>)grad<£> that should model surface tension. However, if the initial condition 
for ip involves an interface of finite width (or, for practical purposes, in any numerical 
scheme), there is no mechanism in floT| ) to avoid spreading of the layer. Hence frequent 
renormalizations of the layer would be required in an actual computation. 



V. A NUMERICAL EXAMPLE: COARSENING 

As a numerical application of the theory, we consider the kinetics of coarsening following 
a quench into the coexistence region of a binary fluid in two space-dimensions (see ref. 
||T5| for a review of coarsening). Numerical solutions of the Cahn-Hilliard equation without 



hydrodynamic effects have been obtained by a number of authors, for both critical [^] and 
off-critical []T7| quenches, while hydrodynamic effects on spinodal decomposition have been 
discussed in ref. [|1T| . We consider here the off-critical region of the phase diagram in which 
the decay of the original state involves the formation of localized droplets. 
Our starting point is Eqs. (|5l|) and fl52[) with, 

/M = -^ 2 + ^ 4 , (62) 

as the homogeneous part of the free energy density, with r and u positive constants. For 
incompressible flows in two dimensions, it is convenient to introduce a stream function (, 
such that, v = ||i — where i and j, respectively, are the unit vectors in the x and 



y directions. We next scale the order parameter by yr/u, lengths by the the mean field 
correlation length £ = Ja/r, and time by the diffusion time r = a/mr 2 . In what follows, all 



variables are assumed to be dimensionless. Taking the curl of ([51]), and restricting attention 



to small Reynolds number (Re = mr/v), so that the inertial term in ([51]) is negligible, we 
obtain, 

A 2 C + C [grad (Aip) x grad(y?)] = 0, (63) 
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where C = 3a^/2r]mr, (77 = pv) plays the role of a capillary number, a = 2ar/3£w is the 
mean field value of the surface tension in this model. The fluid is enclosed in a square cavity 
of side L, with the fluid velocity satisfying no-slip boundary conditions on the walls of the 
cavity. Hence, the boundary conditions appropriate for ( |63| ) are that both ( and its normal 
derivative d(/dn vanish on the boundary. 

The dimensionless Cahn-Hilliard equation reads, 

d t tp + v • grad<£ = —A ip — cp 3 + Acp , (64) 

and is supplemented by the boundary conditions dip/dn = dfi/dn = where p = —[<p — 
(p 3 + Atp). 

The equations (|63|) and (^4|) are solved numerically in a square, two dimensional grid with 
a backward implicit method, and a fast biharmonic solver (further details are given in [|i~8|j). 



We take C = 10, and use the same parameters for the simulation as those of |L7| for a solution 
of the Cahn-Hilliard equation at an off-critical value of the order parameter: L = 256, a grid 
spacing 5x = 1.0, and an initial configuration in which ip is random, gaussianly distributed 
with mean (ip) = l/\/3 and variance unity. We have chosen a variable integration step 
St = 0.01 for < t < 4.0, 5t = 0.1 for 4.0 < t < 664 and St = 0.5 for 664 < t < 3700. 

Our results are summarized in Figures |l] - At short times, drops of the minority 
phase are seen to nucleate, grow, and coarsen so as to decrease the free energy of the 
system. Figure |l| shows the distribution of phases -as defined by the dimensionless order 
parameter- as a function of time, while Figure § displays the corresponding distribution 
for the stream function. The main qualitative difference between our results and those 
for the uncoupled Cahn-Hilliard equation is the flow-induced coalescence of neighboring 
droplets. The times in Figures [I] and |2] are chosen to illustrate coalescence. Figure ^| shows 
large a localized variations of the stream function in the regions in which coalescence is 
occurring. Coalescence is not observed in the purely relaxational Cahn-Hilliard equations. 
There coarsening is dominated by diffusion. 

Figure || graphs the total droplet-perimeter, per unit area, with the droplet boundaries 
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defined as the level set <p = 0. At sufficiently long times growth is faster than the classical 
Lifshitz-Slyozov law (growth law of the form t 1 ' 3 ) for coarsening in the absence of hydrody- 
namic effects. Siggia |nj predicted that hydrodynamic effects would affect coarsening, and 
San Miguel, Grant and Gunton |20| reexamined Siggia's argument for the case of two dimen- 
sions. They argued that the dominant coarsening mechanism in three dimensions (which is 
expected to lead to a growth law in which the average size of the domains is proportional to 
time) is not operative in two dimensions. Droplet coalescence was suggested to be dominant 
instead, and on phenomenological grounds they deduced a t 1 ^ 2 growth law; Figure § is in 
accord with this result. Our results are also in agreement with recent molecular dynamics 



simulations of coarsening in binary fluids ||21|| , but disagree with a previous numerical study 
of the same model in two space dimensions, and inside the unstable region of the phase 



diagram, that obtained a growth exponent of 0.69 [22 
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FIGURES 

FIG. 1. Order parameter values (in grey scale) as a function of time showing nucleation and 
coarsening of the structure. Times have been chosen to highlight coalescence events. 

FIG. 2. Stream function £ in grey scale at the same times shown in Fig. [l]. For the values of 
the parameters used, £ lies approximately between -0.5 and 0.5. The regions of largest variation 
in the stream function are located in the regions where coalescence takes place. 

FIG. 3. Inverse perimeter-density as a function of time for C = 10. The three solid lines 
indicate power laws with the value of the exponent shown. 
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Figure 1, Gurtin et al. 
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Figure 2, Gurtin et al. 
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Figure 3, Gurtin et al. 



